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1.  Introduction. 


The  allocation  of  a -specific  weapon  system  type  to  an  acquire  target  is 
an  important  tactical  decision  in  the  fire-support  process.^  Accordingly,  the 
determination  of  optimal  (or  even  good)  fire-distribution  strategies  for  support- 
ing  weapon  systems  is  a major  problem  of  military  operations  research.  The 
problem  is  of  interest  to  the  military  tactician  so  that  he  may  have,  a clearer 
understanding  of  the  circumstances  under  which  a supporting  weapon  system  (such 
as  artillery)  should  engage  the  enemy's  primary  weapon  system  (i.e.  infantry) 
and  when  it  should  engage  the  enemy's  supporting  weapon  systems. 

Xn  this  paper  we  will  examine  the  dependence  of  optimal  time-sequential 
fire-support  strategies  on  the  form  of  the  combat  attrition  modex.  Previous 
work  by  Weiss  [38]  and  Kawara  [22]  suggests  that  an  optimal  fire-support  strategy 
consists  in  always  concentrating  all  fire  on  one  enemy  target  type  (although 
this  target  type  may  change  over  time) . We  will  consider  a differential  game 
with  slightly  different  combat  dynamics  than  the  fire-support  differential  game 
recently  considered  by  Kawara  [22]  and  show  that  optimal  fire-support  strate- 
gies quite  different  in  structure  than  those  obtained  by  him  may  arise.  More- 
over, the  solution  to  the  problem  which  we  consider  in  this  paper  involves  a 
solution  phenomenon  not  previously  Encountered  in  a Lanchester-type  differential 

*t**fi* 

game:  the  dual  (or  adjoint)  variables  may  be  discontinuous  across  a manifold 

of  discontinuity  for  both  players'  strategies. 

Fire-support  operations  (as  are  any  combat  operations)  are  a complex 
random  process  (see  [26]).  We  will  nevertheless  consider  a simplified  deterministic 

*See  pp.  1-33  to  1-43  oi  [26]  for  a further  discussion. 

See  [38]  for  a brief  discussion  of  the  distinction  between  a "primary"  weapon 
system  (or  infantry)  and  a "supporting"  weapon  system. 

4.  + + 

We  refer  to  a differential  game  as  being  a Lanchestei-typo  differential  game 
when  the  system  dynamics  are  described  by  Luncbester-type  equations  of  warfare 
(see  [34]).  | 
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Lanchester-type  model  in  order  to  develop  insights  into  the  structure  of  optimal 
time-sequential  fire-support  strategies.  H.  K.  Weiss  [38]  has  emphasized  that 
such  a model  of  an  idealized  combat  situation  is  particularly  valuable  when  it 
leads  to  a clearer  understanding  of  significant  relationships  which  would  tend 
to  be  obscured  in  a more  complex  model. 

The  problem  of  determining  an  appropriate  mixture  of  tactical  and  strategic 
forces  (another  aspect  of  the  optimal  fire-support  strategy  problem)  v/as  exten- 
sively debated  by  experts  during  World  War  II.  Some  analysis  details  may  be 
found  in  the  classic  book  by  Morse  and  Kimball  (see  pp.  73-77  of  [27]).  The 
problem  was  studied  at  RAND  in  the  late  1940's  and  early  1950's  (see  [16])  and 
elsewhere  (see  [1]).  It  would  probably  not  be  too  far-fetched  to  claim  that 
this  problem  stimulated  early  research  on  both  dynamic  programming  (see  [2])  and 
also  differential  games  (see  [16],  [20]).  Today  the  problem  of  determining 
optimal  air-war  strategies  is  being  extensively  studied  by  a number  of  organiza- 
tions (cee,  for  example,  [17],  [25],  [29],  [36]).  An  idealized  version  of  A.  Mengel's 
problem  (see  [16])  appears  in  Isaacs'  book  as  the  "War  of  Attrition  and  Attack" 

(see  pp.  96-105  of  [21]).  Discrete-time  versions  of  this  problem  of  determining 
optimal  "air-war"  strategics  have  been  considered  by  a number  of  workers  as 
time-sequential  combat  garnet  [5],  [6],  [15]  (see  also  [7],  [13]).  A related 
problem  has  been  considered  by  Weiss  [.<8]  (see  also  [37]),  who  studied  the  optimal 
selection  of  targets  for  engagement  by  a supporting  weapon  system.  More  recently, 
Kawara  [22]  has  studied  optimal  time-sequential  strategies  for  supporting  weapon 
systems  in  an  attack  scenario  version  of  Weiss'  problem.  Other  recent  work  has 
considered  various  conceptual  and  computational,  aspects  of  time-sequential 
combat  games  [28],  [29],  [30]. 

+ 

See  [33],  however,  lor  a just i t i cat  ton  o'  t Ik*  optimality  < 1 strategies  given 

by  Weiss  [38].  A general  solution  algorithm  is  also  presented  in  this  paper 

[33]. 
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Since  our  work  here  may  be  considered  to  be  an  elaboration  upon  and 
extension  of  Kawara's  fire-support  differential  game  {22],  we  will  review 
his  main  results  and  relate  our  work  here  to  them.  Kawara  [22]  considers  com- 
bat between  two  heterogeneous  forces,  each  composed  of  infantry  (the  primary 
weapon  system)  and  artillery  (the  supporting  weapon  f.ystem).  The  time-sequential 
decision  problem  is  to  determine  each  side’s  optimal  strategy  for  distributing 
its  supporting  weapon  system's  fire  over  enemy  target  types  according  to  the 
criterion  of  the  infantry  force  ratio  at  the  prescribed-duration  attack's  end. 

f 

Kawara  concludes  that  each  side's  optimal  strategy  is  to  always  concentrate  all 
supporting  fire  on  the  enemy's  supporting  weapon  system  (counter-battery  fire) 
during  the  early  stages  of  battle  (if  the  total  prescribed  length  of  battle  is 
long  enough)  and  then  later  to  switch  to  concentration  of  all  fire  on  the  enemy's 
infantry.  He  states  that  this  switching  time  "does  not  depend  on  the  current 
strength  of  either  side  but  only  on  the  effectivenesses  of  both  sides'  supporting 
units"  (see  p,  951  of  [22]).  Moreover,  an  optimal  strategy  has  the  property  of 
always  requiring  concentration  of  all  supporting  fire  on  enemy  infantry  during 
the  final  stages  of  battle. 

Thus,  Kawara  concludes  that  for  his  model  the  optimal  fire-support 
strategies  do  not  depend  on  force  levels.  However,  this  is  only  true  provided 
that  the  appropriate  side's  (in  Kawara's  numerical  example,  the  defender) 
supporting  weapon  system  is  not  reduced  to  a zero  force  level  before  a critical 
time.  Let  us  assume,  therefore,  that  neither  side's  supporting  weapon  system 

1* 

Kawara  does  not  determine  the  optimality  of  extremal  strategies  determined  for 
his  problem  (i.e.  show  that  sufficient  conditions  of  optimality  are  satisfied 
(see  [A])).  We  use  the  work  extremal  to  denote  a trajectory  on  which  the 
necessary  conditions  are  satisfied. 

++  it 

See  the  expression  for  T,.  on  p.  949  of  (22]  and  its  plot  in  Figure  4 of 

[22].  1 
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can  be  reduced  to  a zero  force  level.  For  this  condition  the  optimal  fire- 
support  strategies  are  force-level  independent  and  may  be  expressed  solely  in 
terms  vf  "tiue-to-go"  in  the  prescribed  duration  battle.  The  purpose  of  this 
paper  is  to  show  that  a tactically  realistic  variation  in  the  attrition  equations 
leads  to  a problem  with  force-level  dependent  optimal  fire-support  strategies. 

This  result  has  an  important  implication  for  tactical  decision  making:  optimal 

time-sequential  allocation  of  fire-support  resources  depends  on  not  only  initial 
intelligence  estimates  but  also  on  a continuous  monitoring  of  the  evolution  of 
the  course  of  combat. 

Thus,  the  purpose  of  this  paper  is  to  illustrate  the  dependence  of  optimal 
fire-support  strategies  on  the  nature  of  Lanchester-type  combat  attrition  aqua- 
tions (see  [3*v]).  We  consider  a slight  variation  in  Kawara's  problem  (i.e. 
different  combat  dynamics)  for  which  the  structure  of  the  optimal  strategy  of 
one  of  the  combatants  is  significantly  different  than  that  in  the  original  prob- 
lem [22]:  the  optimal  strategy  of  one  combatant  depends  directly  upon  the  enemy's 
force  levels  and  is  no  longer  to  always  concentrate  all  fire  on  either  the  enemy's 
primary  or  supporting  weapon  system.  Furthermore,  wa  will  show  that  an  optimal 
strategy  in  which  a side  divides  the  fire  of  its  supporting  weapon  system  between 
the  enemy's  primary  (infantry)  and  supporting  systems  can  only  occur  when  the 
enemy's  infantry  has  some  fire  effectiveness  (in  the  sense  of  a non-zero  Lan- 
chester  attrition-rate  coefficient)  against  his  infantry.  The  optimal  strategy 
of  one  side  to  sometimes  split  its  fire  is  very  similar  to  that  which  occurs  in 
a one-sided  (optimal  control)  problem  previously  considered  by  us  [31 . (see  also 
[32])  for  the  optimal  distribution  of  fire  by  a homogeneous  force  in  combat 
against  homogeneous  forces.  In  [31]  the  enemy  consisted  of  two  weapon  system 

^Initial  force  levels  and  the  known  length  of  battle  may  be  sufficient  to 
guarantee  this  for  a given  set  (or  range  of  values  of)  Lanchester  attrition- 
rate  coefficients. 
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types,  each  of  which  undergoes  attrition  at  a rate  proportional  to  the  product 
of  the  numbers  of  firers  and  targets  (referred  to,  for  convenience,  as  "linear- 
law”  attrition).  In  fact,  this  previous  work  of  ours  [31]  was  the  motivation 
for  our  examination  here  of  other  attrition  structures  in  Kawara's  problem. 

2.  Kawara’s  Fire-Support  Differential  Game. 

Since  Kawara’s  fire-support  differentia]  game  is  the  point  of  departure 
for  this  paper,  we  will  review  the  development  of  his  model.  The  reader  will 
find  it  convenient  to  compare  the  mathematical  statement  of  Kawara’s  problem 
(1)  with  the  fire-support  differential  gams  studied  in  this  paper  (2)  in  order 
to  understand  the  dependence  of  optimal  fire-support  strategies  on  the  mathe- 
matical form  of  the  attrition  equations. 

Kawara  [22]  considers  the  attack  of  heterogeneous  X forces  against  the 
static  defense  of  heterogeneous  Y forces.  Both  the  X and  Y forces  are 
composed  of  two  types  of  units:  primary  units  (or  infantry)  and  fire  support 

units  (or  artillery)  . The  X infantry  (denoted  as  X^)  launches  an  attack 
against  the  Y infantry  (denoted  as  Y,) . We  consider  that  phase  of  the  attack 

A. 

which  may  be  called  the  "approach  to  contact."  This  is  the  time  from  the  initia- 
tion of  the  advance  of  the  X^  forces  towards  the  Y^  defensive  position  until 
the  forces  actually  make  contact  with  the  enemy  infantry.  It  is  assumed 

that  this  time  is  fixed  and  known  to  both  sides  and  that  infantry  fire  has 
negligible  effectiveness  against  the  enemy's  infantry  during  this  time.  During 
this  time  the  fire  support  units  remain  stationary,  and  each  unit  has  the  capa- 
bility to  deliver  either  "point-fire"  counterbattery  fire  against  enemy  artillery 
or  "area  fire"  against  the  enemy’s  infantry. 

It  is  the  objective  of  each  side  to  attain  the  most  favorable  infantry 
force  ratio1  possible  at  the  end  of  the  "approach  to  contact"  at  which  Lime  the 
4* 

See  [35]  for  some  insights  into  the  dynamics  of  combat  from  considering  the 
force  ratio. 
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force  separation  between  the  opposing  infantries  is  zero  and  artillery  fires 
must  be  lif.  >d  from  the  enemy's  infantry  in  order  not  to  also  kill  friendly 
forces.  Thus,  the  decision  problem  facing  each  commander  is  to  determine  the 
"best"  distribution  of  artillery  fire  over  time  between  enemy  infantry  and  enemy 
artil.’ery  in  order  to  maximize  the  quotient  of  friendly  infantry  (numerical) 
strength  divided  by  enemy  infantry  strength  at  the  end  of  the  approach  to  con- 
tact. This  situation  is  shown  diagrammatically  in  Figure  1.  The  reader  is 
referred  to  Kawara's  paper  [22]  for  further  details  of  the  model's  development. 
It  should  be  pointed  out  that  this  model  also  applies  to  the  case  of  an  amphib- 
ious landing  and  the  determination  of  the  optimal  time-sequential  allocation  of 
the  supporting  fires  of  Naval  ship  guns. 

Mathematically,  the  problem  may  be  stated  as  the  following. 


x.(t,) 


maximize  minimize 
U V 


/YufM 

\y2(t £>f  * 


with  stopping  rule:  t^  - T - 0, 

dx. 


subject  to: 
(battle  dynamics) 


dt  = 'VVly2' 

dx2 

"dt"  = -(1"v)a2y2’ 

dyl  . 

"dT  = "ublylX2  ’ 

dy  2 

It  = -(1"u)b2X2' 


with  initial  conditions 
o 


i't(t=0)  = and  yi(t=0)  = y±  for  i = 1,2, 


(1) 


*t* 

We  use  capital  letters  to  denote  the  closed- loop  (or  feedback)  strategies  (see 
[19])  of  the  players  and  the  corresponding  lower  case  letters  to  denote  the 
corresponding  strategic  variables  (see  [4]).  A strategic  variable  is  the 
realization  (or  outcome)  of  a strategy.  Thus,  u(t)  = ll(t,£,£)  and  v(t)  = 
v(t,X,£). 


u 


Lgure  1.  Diagram  of  Kawara's  Fire-Suppor 


h-M<- Bgg Wk'HJy’U^i • 4IH v ^ - ‘ ? g? » VgR&mZEgzS' ■w^r; 


and 


x1,x2,y1,y2  ^ 0 (State  Var^^j^^^quaiity  Constraints), 

0 £ u,v  sC  1 (Strategic  Variable  Inequality  Constraints), 
where 

x^t)  is  the  number  of  X infantry  (i.e.  X^)  at  time  t, 
x2(t)  is  the  number  of  X artillery  (i.e.  X2)  at  time  t, 
similarly  for  y^(t)  and  y2(t), 
a is  a constant  (Lanchester)  attrition-rate  coefficient^ 

i 

(reflecting  the  effectiveness  of  Y2  fire  against  X^) , 


similarly  for  b^. 


and 


u(v)  is  the  fraction  of  X(Y)  artillery  fire  directed  at  opposing 
infantry  force-*. 

We  observe  that  for  T<-H»  it  follows  from  the  battle  dynamics  (1)  that 
x1(t)  ,y^(t)  >0  Vte  [0,T].  Thus,  the  only  state  variable  inequality  con- 
straints (SVIC's)  that  must  be  considered  are  x^,y2  ^ 0. 

Kawara's  results  and  conclusions  [22]  have  been  discussed  in  Section  1 

above . 


3.  Another  Model  for  Optimal  Fire-Support  Allocation. 

In  this  paper  we  will  study  a Vitiation  of  Kawara's  [22]  fire-support 
differential  game  (1)  just  given.  We  will  see  that  for  this  problem  the  struc- 
ture of  che  optimal  fire-support  strategy  for  the  attacker  has  a fundamentally 
different  form  than  that  for  (1):  the  attacker  rust  sometimes  split  his  fire 

between  the  defender's  primary  and  supporting  units  in  order  to  "avoid  overkill." 

+ 

See  [10]  (also  [8],  [9])  for  methodology  for  the  prediction  of  such  coefficients 
from  weapon  system  performance  data. 


8 


2£asK-.-„" 


Furthermore,  the  nature  of  this  split  in  an  optimal  strategy  depends  on  the 
allocation  of  the  defender's  supporting  fires. 

Let  us  again  consider  the  attack  of  heterogeneous  X forces  against  the 
static  defense  of  heterogeneous  Y forces.  Each  side  is  composed  of  primary 
units  (or  infantry)  and  fire  support  units  (or  artillery) . The  X infantry 
(denoted  as  X^)  launches  an  attack  against  the  position  held  by  the  Y infantry 
(denoted  as  Y^).  Again,  we  will  consider  only  the  "approach  to  contact"  phase 
of  the  battle.  This  pahse  is  the  time  from  the  initiation  of  the  advance  of 
the  X^  forces  towards  the  Y^  defensive  position  until  the  X^  forces  actually 
make  contact  with  the  enemy  infantry.  It  is  assumed  that  this  time  is  fixed 
and  known  to  both  sides. 

The  X^  forces  begin  their  advance  against  the  Y^  forces  from  a dis- 
tance and  move  towards  the  Y^  position  using  "cover  and  concealment."  The 
objective  of  the  X^  forces  during  the  "approach  to  contact"  is  to  close  with 
the  enemy  position  as  rapidly  as  possible.  Accordingly,  small  arms  fire  by  the 
X^  forces  is  held  at  a minimum  or  firing  is  done  "on  the  move"  to  facilitate 
their  rapid  movement.  It  is  not  unreasonable,  therefore,  to  assume  that  the 
effectiveness  of  X^'s  fire  "on  the  move"  is  negligible  against  Y^.  We  assume, 
however,  that  the  defensive  Y^  fire  causes  attrition  to  the  advancing  X^ 
forces  at  a rate  proportional  to  the  product  of  the  numbers  of  firers  and  tar- 
gets. Two  possible  justifications  of  this  are  as  follows:  because  of  the 

movement  (and  intermittent  concealment)  of  the  forces  and  the  distance 

involved,  the  Y^  defenders  either  (1)  fire  into  a constant  (but  moving)  area 
without  precise  knowledge  of  the  consequences  of  their  fire  or  (2)  when  they  do 
aim  fire  at  X^  targets,  the  time  to  acq  lire  such  a target  is  inversely  propor- 
tional to  the  density  of  forces  and  much  greater  than  the  time  to  kill  an 
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acquired  target.  Under  each  of  these  sets  of  circumstances  the  assumed  form 
of  attrition  has  been  hypotehsized  to  occur  [11],  [37]. 

During  the  "approach  to  contact,"  the  fire-support  units  remain  stationary. 
Each  unit  has  the  capability  to  deliver  counterbattery  fire  against  enemy  artil- 
lery or  "area  fire"  against  the  enemy's  infantry.  In  other  words,  we  assume  that 
each  side's  fire  support  units  fire  into  the  (constant)  area  containing  the 
enemy's  infantry  without  feedback  as  to  the  destructiveness  of  this  fire.  On 
the  other  hand,  the  effectiveness  of  counterbattery  fire  is  not  symmetric  with 
respect  to  the  two  combatants.  We  assume  that  the  defender  has  the  capability 
(for  example,  through  the  use  oi  aerial  observers)  to  sense  when  an  enemy  sup- 


porting unit  has  been  destroyed  so  that  fire  may  be  immediately  shifted  to  a new 

*f*  + + 

target  and  that  fire  is  uniformly  distributed  over  the  survivors.  The  attacker 

however,  either  (1)  does  not  have  the  capability  to  sense  destruction  of  enemy 


> 


c 

f fire  support  units  accurately  (and  hence  distributes  his  fire  uniformly  over  the 


ft 


(constant)  area  occupied  by  the  defender's  fire  support  units)  or  (2)  if  he  does 
have  adequate  fire  assessment  capability,  then  target  acquisition  times  (which 
are  inversely  proportional  to  the  density  of  the  enemy's  fire  support  units)  are 
much  larger  than  the  time  to  destroy  an  acquired  target.  This  leads  to  a 
attrition  rate  proportional  to  the  product  of  the  numbers  of  X2  firers  and 
targets  [11],  [37]. 


Alternatively,  we  may  think  that  the  attacker  has  massed  so  much  artillery  that 
X2  targets  are  always  easily  acquired  by  Y2  once  an  X2  unit  has  been 
destroyed.  Moreover,  it  will  be  assumed  below  that  the  initial  X2  force 
level  is  sufficiently  large  to  guarantee  that  it  is  never  driven  to  zero. 

+ f 

This  assumption  is  not  essential  for  the  structure  of  X2's  optimal  fire  sup- 
port strategy.  A similar  structural  result  may  be  obtained  when  X2'«  attrition 
is  the  same  form  as  that  for  Y2.  We  have  made  the  above  assumption,  moreover, 
so  that  the  resultant  attrition  model  is  most  similar  to  Kawara's  [22]  but  yet 
yields  significantly  different  results  for  the  attacker's  fire  support  strategy. 
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It  is  the  objective  of  each  side  to  attain  the  roost  favorable  infantry 
force  ratio  possible  at  the  end  of  the  "approach  to  contact"  at  which  time  the 
force  separation  between  the  opposing  infantries  is  zero  and  artillery  fires 
must  be  lifted  from  the  enemy's  infantry's  position  in  order  not  to  also  kill 
friendly  forces.  Thus,  the  decision  problem  facing  each  side  is  to  determine  the 
"best"  distribution  of  artillery  fire  between  enemy  infantry  and  artillery  over 
time  in  order  to  maximize  the  infantry  force  ratio  at  the  time  of  contact  between 
the  two  infantry  forces.  This  situation  is  shown  diigrammatically  in  Figure  2. 

The  above  assumptions  lead  to  the  following  differential  game  with  an 
attrition  structure  slightly  different  than  that  in  Kawara's  problem  [22]. 

xi(cf>i 


maximize  minimize 
U V 


\y1(tf)  f 


with  stopping  rule:  t^  - T = 0, 


subject  to: 
(battle  dynamics) 


dx1 

dt 

-ailXlyl  - ' 

dx_ 

4 

dt 

-(l-v)a2y2. 

dyl 

dt”  = 

-uViv 

dy2 

£.  __ 
dt 

-(l-u)b2y2x 

Va12Xly2' 


(2) 


with  initial  conditions 


x^(t=0)  = x°  and  y^(t=0)  = y^  for  i = 1,2, 


and 


Xl,X2,yl,y2  ^ ® (State  Variable  Inequality  Constraints), 

0 * u,v  ^ l (Strategic  Variable  Inequality  Constraints), 
where  all  symbols  are  (essentially)  Liu:  same  as  defined  above  for  problem  (1). 
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We  observe  that  for  T < +•  it  follows  from  the.  battle  dynamics  (1) 
that  x^(t),  y^(t),  and  y2(t)  >0  Vt  c [0,Tj.  Thus  the  only  SVIC  that  must  be 
considered  is  X2  5*  0.  However,  let  us  assume  that  the  force  level  of  the  attacker's 
artillery  is  never  reduced  to  zero.  In  other  words,  we  consider  the  special  case 
in  which  x°  and  T are  such  that  X2<T)  > 0. 

4,  Characterization  of  Optimal  Fire-Distribution  Strategies  for  the  Supporting 
Weapon  Systems, 

It  snould  be  clear  that  in  (2)  above  we  have  a;Q»ay2’a2*bl*b2  > Although 

the  results  of  A.  Friedman  [14]  concerning  existence  of  value  do  not  apply  directly 

to  our  fire-support  differential  game  (2) , they  do  apply  to  a suitably  modified 

dx2 

version.  If  we  were  to  consider  a version  of  this  problem  with  -j^—  * -(l-v)a2Y2+ 
where  ^ > 0,  then  it  may  be  shc*n  (see  pp.  210-230  of  [14])  that  this  "modified" 
differential  game  has  value  and  that,  a saddle  point  exists  in  pure  strategies  (see 
pp.  2%-235  of  [14]).  We  will  now  develop  the  basic  necessary  conditions  of 
optimality  for  (2). 

For  xi»x2»yi»y2  > t*ie  Hamiltonian  for  (2)  is  given  by  [12] 
H<t>&>X.>H>aaU,v)  = -Pi(a11Vl+va12Xly2)  ' P2a2(1_v)y2 

" VblylX2  ' q2(1_u)b2y2X2’  (3) 

where  we  have  adopted  the  following  correspondence  between  state  and  adjoint 
variables : 

state  variable  dual  variable 


for  i = 1,2. 


The  adjoint  system  of  differential  equations  for  the  dual  variables  is 


ipi 

9H 

with 

dt 

’ aiiyipi  + 

V*a12y2Pl 

dp2 

3H 

dt 

" *2 

” “*biyiqi 

t (i-u*)b2y2q2 

with 

dql 

d! 

3H 
" 3yl 

= anxipi  + 

u*bix?ql 

with 

dq2 
d~  = 

3H 

" 3y2 

= v*a12XlPl 

+ (l-v*)a2P2  + 

(l-u*)b2x2q 

<^>2 


The  results  of  Berkovi'.z  [3]  say  that  H,  £(t),  and  <^(t)  are  continuous 
functions  of  time  except  possibly  at  manifolds  of  discontinuity  of  both  U*  and 
V*  (see  Section  4.3  below). 

When  x^,x2,y^,y2  > 0,  the  extremal  strategic-variable  pair,  denoted 
as  (u*,v*),  is  determined  by  the  max-min  principle.  Hence,  we  consider 


maximize  minimize  H(t.,x^,{>,c^,u,v)  , 
0 — u— 1 0*v£] 


so  that 


u*(t)  = 


1 for  Su(t)  > 0, 


0 for  S^(t)  < 0, 


where  the  ’J-switching  function,  S^Ct),  is  given  by 


Su(0  = bl("ql)yl  ~ b2('q2)y2> 


v*(0  = 


1 for  S (£)  > 0, 


0 for  SvU)  < 0, 
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. a»U^l*'!lvr'”  'y^ss^jpsri^sgptKSff^'^fv; 


where  the  V-switching  function,  Sv(t),  is- given  by 


Sv(t)  = allplxl  ' a2P2‘ 


It  is  readily  shown  that 


P1(t)x1(t)  = constant  = p^(T)x. (T)  = , 

3 y, 


dt  *qlV  = aU  f yl^  > °» 


f " "a2(1_u*)5M(t)  " a2biqiyr 


We  must  further  investigate  the  possibility  of  singular  subarcs  (see  [31] 


or  Chapter  8 of  [12]).  Let  us  furst  show  that  it  is  impossible  to  have  a 


V-singular  subarc.  In  other  words,  v*(t)  must  be  0 or  1 almost  everywhere 


in  time.  The  impossibility  of  a V-singular  subarc  is  established  by  showing 
ds 


that  > 0 for  all  te  [0,T].  It  is  clear  that 


(l-u*)Su(t)  & 0 for  all  t £ [0,T]. 


Considering  (13)  and  the  fact  that  q^Tjy^T)  j < 0,  we  see  that 

yf 

ql(t)y^(t)  < 0 for  all  t £ [0,T],  whence  follows  the  assertion  via  (14). 


It  is  possible,  however,  to  have  a U-singular  subarc  on  which  -j—  -•  0 


(or,  equivalently,  ( t ) = 0)  for  a finite  interval  of  time.  There  are  two 


cases  to  be  considered:  (1)  v*  = 1 and  (2)  v*  = 0. 


4.1.  U-Singular  Subarc  on  Which  V*  = _1 . 


When  v*  ° 1,  it  is  readily  computed  that 


, r 

P‘11,  . , , 

r ff  (aiil,iyrai21,2y2)- 

V\i 
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gr-y» y» ■n^jSjgg^ .j—  T ,?.y- 


i 


and 


d^S 


dt 


7~ 


X5 i ^a]  1blyl)'l*bl  " (ai2b2y2^1"U*^b2^’ 


(17) 


3H 


Considering  (9),  the  requirement  that  ~-  = 0 yields  the  first  condition  for 


OU 


a U-singular  subarc  with  V*  = 1 


blVl  ■ b2V2- 


(18) 


Considering  (16)  and  (.18),  the  requirement  that  —-(•—)  --  0 on  a singular 


8H 


subarc  on  which  — = 0 for  a finite  interval  of  time  yields  the  second  condi- 
tion for  a U-singular  subarc  with  V*  = 1 


M 

l 

I 

i 

l 


% 


I 


| 

V 

£ 

& 


allVl  ■ a12b2y2  OS) 


On  a subarc  on  which  the  first  and  second  conditions  for  a singular  subarc 

2 

hold  we  additionally  require  that  —7  (~)  = 0 so  that  (17)  yields  the 
singular  strategic  variable  value  required  to  keep  the  system  on  the  singular 
sub arc 


u*(t)  = 


bl+b2 


(20) 


Checking  the  generalized  Legendre-Clebsch  condition^  [23],  [24]  ^{--7^--]}  ^ 0, 

duldt  l9uy J 

we  find  that  on  a subarc  on  which  (18)  and  (19)  hold  we  have 


f 

rx,'i 


&E7©*  = H}(x2)?'Wl(bl)2yl  + a12(b2)2y2}  " °* 


£ 

'H 

I 


We  may  write  the  equation  of  the  U-singular  surface  (see  p.  683  of  [31]  as 


yi  a12b2  f 

__  = for  v*  = 1. 

y2  allbl 


(21) 


I’ll  is  is  a necessary  condition  for  optimality, 
equivalent  condition  (seo  [18]). 


R.  Isaacs  [31 | gives  an 


4.2.  U-Singular  Subarc  on  Which  V*  - 0. 

When  v*  = 0,  it  is  readily  computed  that 
dS  xf 

dt“  = ” “f  (allbiyra12b2y2)  ” b2y2Sv<t)’ 

yl 


and 


(22) 


d2S  dS 

dr^  ’ •“*blx2  dT  + a2b2y2(-u*Su(t)+b2q2y2+p2*2[u*b1-a-u*)b2)}, 


(23) 


so  that  the  first  and  second  conditions  for  a U-singular  subarc  with  V*  *=  0 
are,  respectively,  (18)  and 


allVl  = a12b2y2  + b2y2 


f <“Sv<t)}* 


(24) 


It  should  be  noted  (see  [lb])  that  the  above  singular  surface  exists  in 
£ - ^ space.  It  is  convenient  to  write 

yl  a12b2  . b2y)  r „ , ,, 

— = ~L-X~  + (t)}  for  v*  = 0. 

y.  a,.b.  , f v 

2 11  1 anVl 

The  singular  strategic  variable  value  is  given  by 


(?5) 


u*(t)  = 


b2  I 

fl  V21 

lbl+b2J 

W 

(26) 


The  requirement  that  u*  £ 1 yields  that  on  a b-singular  subarc  with  V*  * 0 
we  must  have 

b2(“q2)y2  * blP2X2*  (27) 

It  is  readily  checked  that  the  generalized  Legendre-Clcbsch  condition  is 
satisfied. 

4.3.  Discontinuity  of  Adjoi nl  Variables  Across  Manifold _of  Discontinuity  of 
Both  U*  and  V*. 

It  is  convenient  to  introduce  the  hackwjrds  time  t defined  by 


T = T - t. 


(28) 
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From  (20)  and  (26),  we  see  that  u*(t)  must  change,  in  general,  discontin- 
uously  from  b2^bl+b2^  to  b2^bi+b2^  whenever  v*(t) 

changes  from  1 to  0.  Let  us  consider  the  totality  of  trajectories  on  which 
this  happens.  The  locus  of  points  in  the  t,x^-  space  for  such  simultaneous 
switches  is  then  a manifold  of  discontinuity  of  both  U*  and  V*.  Across  such 
a manifold  the  adjoint  variables  need  not  be  continuous  (see  [3]). 


J,et  *»  tv(j<,£)  denote  the  backwards  time  at  which  v*(t)  changes 


from  1 to  0.  For  future  purposes,  it  will  be  convenient  to  consider  a 


simultaneous  switcli  with  u*  changing  from  the  singular  control  b2^bi+b2^ 


to  1.  Then  the  manifold  of  discontinuity  of  both  U*  and  V*  is  given  by 

F(t,x,£)  = t - T + tv(x,£)  = 0, 


GW  = auVi  - a12b2y2  = °* 

Across  the  manifold  of  discontinuity,  we  have 


T,  +, 

u.  <TV> 


T,  3F  9G 

H (Tv)  “ p " 0 ’ 


T.  +. 

a.  <0 


T,  -v  9F  DG 

a-  <V  - p 3^  - 0 ■ 


"(<) 


9 F . 9 G 

Tv>  + " Tf  + 0 37  • 


T . +.  T. 

’ (\)  " P (tv) 


V 

p ST  ’ 

rv 

(30) 

9x 

a^  + oanbi  • 

(31) 

3t 

Fyj  " oa12b2  • 

(32) 

• P • 

(33) 
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Considering  (9)  and  (11),  it  is  readily  shown  that 


V\>  ' 0(all<bl)?>'l  + a12tb2>V  + p(Vl  3y‘i  ' b2y2  ^ ) ' <34) 


and 


+ 3t  3t 

S,(Tv>  ■ ~P^ai2xl  'taj  ' a2  tej’- 


(35) 


Recalling  that  u*(tv>  = b2^bl+b2^’  u*(t*)  * 1,  v*(tv>  = 1,  and  v*(t*)  = 0, 

3x 

we  may  substitute  (30)  through  (32)  into  (33)  to  obtain  for  a.-x.y  + 

3Ty  3tv  11  1 * 

a2y2  3^  + biyiX2  3^  ^ 1 


P = 


aU(bl)2ylX2° 


3t  3t  3t 

(l-anxiyi  g—  - a2y2  7—  - 


(36) 


Then  we  may  write 


VTv>  ' °{aU<bl)2yl  + a22(b2)2y2 


and 


3t 


3t 


*n<bi)  V2(biyi  3^  - V2  ^ 


3t  3t  3t 

(l-anxiyi  yx~  ~ a2y2  Jx2  ~ bly)X2  3y^' 


I 


S (t+) 

V V 


•all<bl)V2(a12xl 


3t 


3t 


v v. 

3x^  a2  3x2 


3r 


3t  3t 

(1_aUXiyl  3*"  “ a2y2  3x~  “ blylX2  Dy^ 


(37) 


(38) 


5 . Synthesis  of  Kxt vernal  Strategic- Variable  Pa i_r . 

By  the  synthesis  of  the  extremal  strategic  variable  pair  we  mean  the 
explicit  determination  (using  the  basic  necessary  conditions  oi  opt  imali  ty) 
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of  the  time  history  of  the  extremal  strategic  variable  pair  (u*,v*)  from 
initial  to  terminal  time  (see  (21]  and  also  ( 31 ] — [ 33 ] ) - The  basic  idea  is 
to  trace  extremals  backwards  from  the  terminal  manifold  (where  boundary  con- 
ditions for  tiie  adjoint  va"iables  are  known)  in  such  a way  to  guarantee  the 
satisfaction  of  the  initial  conditions.  Thus,  it  is  convenient  to  introduce 

the  backwards  time  t defined  by  (28). 


5.1.  Extremal  Transitions  in  Strategic  Variables. 

It  seems  appropriate  to  examine  what  are  the  possible  transitions  (ot 

changes)  in  each  strategic  variable  as  we  work  backwards  from  the  end  (i.e. 

3S 

as  t increases).  It  has  been  shown  above  that  < 0 for  all  x C [0,T]. 

Considering  the  boundary  conditions  (4)  and  (5)  for  the  adjoint  variables, 

it  follows  that  S (t=0)  > 0.  Thus 
v 


v*(t)  = 


1 for  0 * t ^ t , 
v* 


0 for  x < t. 
v 


It  will  be  convenient  to  refer  to  that  phase  of  the  planning  horizon  during 
which  v*(t)  =0  as  V-Phase  I (i.e.  0 £■  t < T - t ) and  to  that  during 

which  v*(t)  = 1 as  V-Phase  II. 

Extremal  transitions  in  u*  for  increasing  x as  shown  in  Figure  3. 

Thus,  this  figure  shows  what  changes  we  might  expect  to  obs.erve  in  u*  as 

we  follow  an  extremal  backwards  from  the  end  of  the  planning  horizon  at  x = 0. 

dS 

During  V-Phase  II  when  v*  = 1,  2^2y2^  wit^  S (x=0)  = 

f f d fyil  1 

b.x./y.  > 0.  When  u*  = 0,  then  3 <0.  During  V-Phase  I when  v*  = 0, 

* 1 1 f dTly2J 

d$u  /■x^'i 

- — = — r (a. .b.y.-a. „b0y„)  + b^y.S  (x).  During  both  phases,  the  singular 

ax  t nit  1/  t i.  t.zv 

yl 

subarc  may  be  exited  with  either  u*  = 0 or  u*  = 1.  Once  u*  becomes  0, 
it  remains  this  way.  The  above  statements  will  bo  further  justified  below. 


+It  should  bo  kept  in  mind  that,  for  example,  u*(t)  = I'Mt.x^v) 


P2*2 


for 


for 


v*  = 0, 


v*  = 1. 


Figure  3.  Extremal  Transitions  in  u’>  for  Increasing  t 


5.2.  Ex  t re mal  Syn the s is_  for r * __<_ t ^ . 

From  tin*  above  we  have 


so  that  by  (8)  we  have 


S (t=0)  = b 4 > 0, 

>1 


u*(t)  = 1 for  0 £ r £ T , 

u 


where  t is  the  smallest  zero  of  the  equation  S (t=t  ) = 0.  If  the  U- 
singular  subarc  is  reached  in  V-Phase  II  (see  Section  4.1.  above),  then  let 
us  denote  the  backwards  switching  time  at  which  u*  changes  from  1 to 

(bi"*'b 2)  as  t*.  Clearly,  it  is  necessary  that  t*  < for  this  singular 
subarc  to  appear  in  the  solution.  Thus,  in  genera],  there  are  two  cases  to 


be  considered: 


(1)  T*  < x , 
u v 


(2)  T*  T . 

U V 


In  this  paper  we  will  consider  only  the  former  case,  with  the  latter  one  follow- 
ing along  the  same  general  lines  of  development.  We  therefore  assume  that  a-Q» 


, , f f . f 

2*  a2 * x2 * 


f * 

r.  are  such  that  t < t . We  will  give 
-2 u v 


numerical  results  for  this  case  below.  Moreover,  in  all  our  numerical  computa- 
tions we  have;  only  encountered  this  case. 

5.2.1.  Ex i retrials  Near  the  Tormina ] Mani fold . 

dSu  " y} 

Recalling  (16),  we  see  that  — — > 0 (<0)  if  and  only  if  — > (<) 

a12b2  ‘ ' ^ 

r~  . Considering  (40)  it  is  clear  that  S (t)  > 0 for  v*  = 1 when 

•Vi  u f 

a12b2  a12b2 
> . However,  S (r)  may  change  sign  when  — 7 < — ~ . The  U- 

y2  ll  l y 2 ill 


singular  subarc  occurs  when  both  S (t=v->)  *-*  0 and  a..b,y.  - a.„b„y,,  at 

u 11  11  1 1 12  2 2 


T = t*.  Thus,  t*  is  the  smallest  root  of 
u u 


22 


i 


mmw 


— r + T * + 

. f u 
blX2 


1_ 
" f 
'1X2 


b,x,  a y 


f 

11' 1 


~bi  V* 

! 2 u _ 
e * 0. 


(42) 


If  y^  is  given,  then  Su(t=t*)  = ® and  allbl^l  = a12b2y2  be  combined 

to  yield  the  value  of  y^  required  in  order  to  reach  the  U-singular  subarc 

f f f f bi  ^aT  i y f—b i xo ) 

(denoted  as  y *).  Thus,  for  a-.y,  ^ b,x.  we  have  y*  = -—  — — — — ---- — . 

* 11  i ll  2 b0  , r .. 

2 a12(l-blV*) 

(Other  results  are  given  below  in  'fable  I.)  We  denote  the  corresponding  ratio 

f f fyl 

of  y.  to  y * as  — r 

1 2 (y^. 

When  ^U^T=TU^  = 0 with  < ai2b2y2*  ^°H°WS  tbat  Tu 

the  smallest  root  of  the  transcendental  equation 

f 


bl  - 


anyi 


3t‘ 


f u f 

, f ,allyi  blVu  n 
• al2b2y2Tu  + ~T~  e = °- 
X2 

i,  f 


(43) 


It  may  be  shown  that  t— - > 0 where  r = y,/y0.  This  latter  result  is  useful 

dr  12 

In  proving  the  following: 


THEOREM  1:  Assume  that  t > t*.  Then,  u*(t)  = ] on  any  extremal 

v u f f 


1 


as  long  as  v*(t)  = 1 for  — ^ 


yl^ 


PROOF:  The  proof  is  by  contradiction.  Let  r = yjVyjj. 

£ (a)  Assume  that  we  could  have  a switch  in  u*(t)  (with  v*(t)  = 1) 

M ...f. 


for  -j  > 
y2 

with 


^2 


In  other  words,  we  can  find  t such  that  S it*!  ) - 0 

u u u 


S (t)  > 0 for  0 i t < i , 
u u 


(44) 


for  “•  > 


.V 

y2 
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»-  ■ * jV'*-  1 


^ yi 

(b)  Consider  r = — 

y 

3? 


yl' 


2 y2 


+ e with  f.  > 0 and  such  that  t < t . 

u v 


u 

Then  it  may  be  shown  that  — ^ >0.  In  particular  ---a 

dr  3r 


vN*  > 0.  This 


V- 

*_L 


implies,  however,  that  t > t*  for  r = r. 

u u 

u*  » 1 


(c)  Observe  that  J , . for  0 * t f t <t  so  that  y,  /y„  = 

tv*  * 1 uv  12 


f . c biV  „ 

yl'y2  6 ‘ Hence» 


yi  , . *12b2  f 

t2  (t  V " ^ for  r = r’ 


'45) 


dS 


since  then  t (r=r)  > t*.  It  has  been  shown  above  that  — > 0 for  y./y.,  > 
u u ._  dr  12 

do 

a12b2^allbl^  * ^us»  (45)  implies  that  -j^—  (t«t  ) > 0,  and  hence 


0 = S (t=t  ) > S (t)  for  x C (t  -6  ,t  ). 
u u u u u 


(46) 


This  last  statement  (46)  is  a contradiction  to  (44)  , and  the  theorem  is 
proved.  Q.E.D. 


Other  results  are  obtained  in  a similar  fashion. 


5.2.2.  Field  Construction. 

For  a given  set  of  terminal  values  x|,  x^,  y^,  and  y^  an  extremal 
, may  be  traced  backwards  from  the  terminal  manifold  by  a backwards  integration 

i 

| of  the.  state  and  adjoint  equations  combined  with  (8)  and  (10 ) (also  (20 ) or 

(26)).  By  varying  these  terminal  values  the  entire  field  of  extremals  (see 
j p.  128  of  [12])  may  be  obtained. 

; fhe  various  types  of  extremals  that  may  occur  in  the  field  of  extremals 

( are  shown  in  Figure  4.  This  figure  is  representative  of  all  our  numerical 

results  for  t*  < t (see  Section  5.2.3  below).  Pertinent  information  oon- 
! u v 

i cerning  each  type  of  extremal  is  given  in  Table  I. 

I 

t 

! 

i 
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Figure  4.  Identification  of  Various  Types  of  Extremals  for  Which  Information 
is  Given  in  Table  I. 


Table  I.  Extremal  Trajectories  for  Fire  Support 

* 

Problem  with  r < t . 

u v 


u*(t)  = 1 


u pso:  \ . for  °STStu  with  T=  'f 

' v*(t)  =1  y2  'y2J 


f f . 

h fM* 

j = LfJ 


Tu  is  the  smallest  positive  root  of 


1 , * 1 1 } “blx2Tu  _ 

f+Tu  ,f~  fje  " °* 

blx2  blx2  allyl 


with  the  following  bounds  established: 


* f . . f 

£°r  auyx  > b,x2: 


- f . f 

for  a^^b^: 


for  auy[  < b^: 


1 * 1 
f < Tu  < ‘ f ’ 
aliyi  blX2 

T*  . — i — 

Tu  . f * 
blX2 

1 < / < -1 
, f u f 

blX2  Uliyi 


f , , f f*  bl  (ailyl“blX2) 
for  «uyi  * V2:  y2  - , 

2 a]2  1 blX2Tu^ 


r f f f*  allbiyl 

for  a^-b^:  y2  - — e. 


PM  aJ2  a13y] 

S (t)  » a0b.  -4  { 
v 2 1 f ) a b f>? 

yx  ( L2  i (b1x2)  . 


f 

f , f * 

+ [auh.  „ 

,L  . W.L  ebiV 

U / 

r1  *9 

' (b.xh7  i 

w ■ 


Tabic  I.  (cont.)  - 1 


1.  PgyJ  (concluded) 


Let 


Sy(T=Tu).  Also,  on  PgQ  we  have 

f v f 

f f anyi  bix9T 

x1(t)  = Xl  ex5){a12y2T  + (e  -i)}  , 

blx2 

x2(t)  = x2f  , 

. f 

f blV 

y^t)  = yj  e 1 2 , 


and 


y2M  - y2  . 


Px(t) 


P2(t) 


qx(t) 


f . f 

l f anyi  bix?T 

-7  exp  t-[a12y2T  + (e  -1)]}  , 

yl  blx2 


= b. 


V 


vy 


t 

anyr 

{ 11  - — -p 
b!X2 


T + - 


f . f 

>llyl  blx2T 

11  1 (e  1 1 -1)}  , 


(b^)2 


rXl 

Kyl' 


{--u  + 

blX2 


11 


' -blx2T 


x v ^ 
°1X2  yl 


q0(T)  = -a 


12 


rx, 


y 
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— - ~ * /•>>**■  •*  - - 


Tabic  1.  (corn.)  - 2 


n , »*M  - b2/(b1+b2)  * * 

2-  P«-  'uS,i'v 

* II 

where  tu  Is  determined  in  1.  On  Pgj  we  have 


Su(t)  ■=  0, 


allbiyl  * a12b2y2* 


* * 

t is  the  smallest  positive  root  of  S (t*t  ) = 0,  where 
v r v v 


v - < * wrll) 

r . f a..n 


U*fyil  . ^1(T-T*> 


(0x2)2 


exp[0x2(T-Tu) ] | , 


7* 

with  0 « b^b2/(b^+b2) . An  upper  bound  on  tv  is  given  by 


Tv‘  *12/<*2bl)' 


Also,  on  Pgl  we  have 


u*  , ..fir  „ f / *\ 


. - u*  . f-llyl  +a12y2l r e*2<T”Tu)  1 i 
<l(T>  *-  x exp  < j |_e  -1J  J , 

' L 0x- 


x2(t)  ■=  x2  , 

, , u*  0X2(t'tu) 
y2  (t  ) B y2  e 

f 0X2<T-T*) 
y2^xj  “ y2  e 


, , U*  / all>T+e12y2ir°x2(T-V  ,11 

P1(t)  - Px  exp  j -£ 1 [e  -lj( 

1 L ex  j 


with  P;*  - , 

xi  V 


....  u*  u*fXlU  all  [”  U*fyl^  all"|.  * all  6x2fT"1u) 

P2w  - P2  - Vl  ^)!-y  - [-!  [7]  ax?  <~v  - (0xf)2 


* 


f * 


u*  "®x2 (T_TU)  PM*n  , -0X2<T-Tu)l 

i(t)  “ *>1  e - Nrn1  ■ e > • 


yl  0x2 


, . u* 

<)2(t>  “ <l2  e 


-°x2f(T-/)  rx[^ 

" f.  - f I 

V °X2 


.2 

-1-e  ) • 
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Table  I.  (cone.)  - 3 


3 

3*  PA1* 


U*(t)  - 1 
V*(T)  - 1 


for  0 s t t with 


yf  (y\* 

y2  y2 


t ia  the  smallest  positive  root  of  Sv(t*Tv)  * 0,  where  Sv(t) 
la  given  in  1.  An  upper  bound  on  is  given  by 


T “ 


12 


v a^ 


It  haa  been  shown  that  Su(t)  > 0 for  0 * x * t . The  solutions  to  the 
state  and  adjoint  equations  are  the  saae  as  those  for  P**  given  above. 


I / U*<T)  * 1 I 

4.  P7t:  ( for  r 

M W>-0 


We  have  that  S (t)  <0  for  x > x and  that 
v v 


s*  ">  • vk<’> + 


BKVi 


(T)-a12b2y2>- 


Also,  on  we  have 

dx 


1 v 

dT  “ anxiyi  wlth  xi(T*V  ‘ V 


x2<t)  ■*  x2  + a2y2^T“Tv^’ 


a,b.y. 


Yl(T>  « expfbjX^T-^)  + -2~—  (T-Tv)2), 


and 


y2<T>  *y2* 


dp 


~ - -anylPl  with  Pl(-Ty)  « py. 


22 

dt 


~blyiql  With  p2(T"Tv)  " p2’ 


dq 

dt 


~ * _8ii  *T  ’ bix2qi  vith  V^V  • V 


dt”  " -a2p2  With  q2(T“TV>  * V 

We  have  not  been  able  to  develop  solutions  in  terras  of  "elementary"  functions 
to  the  equations  for  x^,  p^,  p^,  , and  q£. 
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Tabic  I.  (cone.)  - 4 


u*(x)  - 1 


>•  # { 

v*(x ) - 1 


for  T_,  £ X * X 
SL  v 


T is  the  smallest  positive  root  of  S (x*=x  ) - 0,  where 
v r v v 

y!  '°1X2'  X1 


. f.  II, 
t (T“TC1 ) 


ill  , iiv  ii  iV  si/l 
“F  (t-tsl)  “ 7 f77  e /• 


bjX^  (b^)2 


Again,  an  upper  bound  on  xv  is  given  by  a^2^a2bl^‘  lc  bas  been  s^own 

that  S (t)  > 0 for  x**  < x * x . Also,  on  P**  we  have 
u SL  v A2 

SLII  . ,n,  .. 

, . SLII  ( SLI1,t  II.  . aliyi  rblX2(  SL}  .] 
x1(t)  - xx  axpja12y2  (x-x^)  + —— — [e  -lj 


x2(x)  «=  x2, 

. . SLII  blX2^T_TSL^  * 
* yl  e 


, , SLII  f SLII _II,  ailVl  f 1 
Pl(x)  « Pl  exPj-a12y2  (x-xSL)  - — j-[e 


SLII  , f,  II, 
anyi  rVz^SL* 


with  p 


SLII  1 1 


SLIllfj 
X1  yl 


, , SLII  , SLII fXll I 3 

p2<x)  - P2  - b y hr  -- 


Rift*  frrf-i 

Vyl  blx2'  L \ 


4 11  / II, 

+ ~f(T‘T  SL} 
b2X2 


. f/  II, 
an  blX2  T-TSL  ) 

. II,  f . f,  II, 

, . SLII  “blX2<T'TSL)  fXl)  311  /,  _blX2(T"TSL)\ 

q1<T)  qi  e " T ~T 

S'  blX2 

f 

/ \ SL1I  fXll  , II, 

q2(T)  - q2  " a12[-FJ(T-Xs1p. 


S. 


•*  tSL 


r-’t''^'  ••  *t  -v» 


t ' v>iC'^5  i',*S’4?1 


u*(“0  * 1 


6*  P*0:  { for  t * t * t 

A2  « V*(T)  - 0 


Table  I.  (cont.)  - 5 
I 


Results  are  similar  to  those  for  above  in  4. 


u*(t)  * 1 
v*(t)  - 0 


1 ( 1 

K ^ for  1 Jfc  T 4 T 


Results  are  similar  to  those  for  P.,  above  in  4. 

Al 


x / u*(t)  * b-/(b.4b  )*(l-q  y /(p  X.))  * u * , 

8.  P*:l  2 12  22  22  S for  t a.  t s t1 

S2  « v*(t)  - 0 v u 


As  usual,  we  have  that  S (x)  < 0 for  x > x . In  order  for  a 
— v v 

* 

U-singular  subarc  to  be  possible  for  x s xy  the  following  condition  must 


hold  at  x - x 


*- 


blp2(Tv")x2<Tv)  51  b2('q2(Tv"))y2(Tv)- 
Also,  on  Pg,  we  have* 


ii. 

dx 

allVl 

with 

Xj^x1 

‘0 

V* 

s y 
X1 

dx2 

dx 

a2y2 

with 

x2(tr 

V* 

a v 
X2 

dyl 

★ 

with 

yi(r 

V* 

dx 

uSblx2yl 

V 

= yl 

dy2  „ 

* 

V* 

dx 

(l-us)b2x2y2 

with 

y2k' 

‘Tv> 

" y2 

and 


dPl 

dx 


*h. 

dx 


A_h 

dx 


“allylPl 

with 

pi(T"0 

V* 

" pi  * 

“blylql 

with 

P2U”0 

V* 

" P2  • 

< *b 

’all  1 _ USblX2qi 
yl 

with 

V* 

“ ql  * 

-a2P2  - (l-u*)b2x2q2 

with 

q2(T=x*) 

V* 

- q2  • 

*A  further  discussion  of  the  continuity  of  the  adjoint  variables  is  to  be  found 
in  Section  5.2.4.  below. 


Table  I.  (cont.)  - 6 


9.  P 


u*(t) 


V*(t) 


for  0 s.  t & x with 


^ < (¥ 

tr  mi  ' 


is  the  smallest  positive  root  of 


k-^)- 

V V 


hih  Vk 


a12b2^2Tu  + * 


f f 

It  should  be  noted  that  > 0,  where  r = y^/y*.  lt  may  be  shown  that  for 


yi  [yiV 
0 ‘ 1 < -7 

v W ^ 


a12y2 


2}  * 
< t < T , 
f u u 


where  the  determination  of  Tu  is  given  in  1.  We.  also  have  that 


Tu(rj)  < Tu^r2^  for  ri  < r2  ^xi  and  x2  held  constant).  The  solutions 


to  the  state  and  adjoint  equations  are  the  same  as  those  for  given  above. 

DU 


Let  Sy(T=Tu)  = S“,  Pl(t=Tu)  = p”,  etc. 


Table  I.  (cont.)  - 7 


TT  / U*(t)  * 0 

10.  Pft9:  { for  T 

82  * V*(x)  - 1 


u 


y 

It  follows  that  for  all  T > Tu  we  have  S t (T)  < 0 and  ~ (t)  < 


u 


ai2b  2^allbl^*  Tv  is  tbe  smaHest  positive  root  of  S (t»t  ) » 0,  where 


Sv(t)  is  given  by 


yl  ih2X2Y 


uPl 


a 


% -f  + 


12 


2lxfJ  ' b xf 
X1  b2X2J 


<T'T„>  -7rVle,'Plb2X2(’'Tu)1}- 

0>2X2)Z 


JI 


Also,  on  Pg2  we  have 


f . i,  m ; 

XjCT)  - xx  exp{a  y"(T-T  ) + [e  2 2 " -1J>, 

b2x2 


x2(t)  - x2, 


y-po  “ yl» 


£ £ 

y2(T>  “ y2  exP^2X2^T_Tu^  * 


and 


Pi  CO  “ P^  exp{-a^yb(T-Tu) [e  2 2 u -1]}  with  p“  = 

b2X2 


f 

1 

rxii 

u 

f 

X1 

V 

P2(T)  - p“  - b2y2 


A] 

/ S12 

>!] 

+ _ln 

>(b2x^ 

2 xfJ 

L x! 

+ . f 
b2X2- 

au  bxh-,) 

(t-Tu>  - nr 0 


(b2x2>2 


qx(T)  - q2  - an 


f 

^x^ 


fxiO  a 


1 0 


fa,. 


. ryfil  -b,4(T-T..h 


■it 

Wr- 


I 


I 


u-  Tlr  { 


u*(t) 

V*(t) 


Table  I.  (cont.)  - 8 


II 


for  x„T  s t <t 
SL  v 


Tv  is  the  smallest  positive  root  of  Sv(x=xv)  * 0,  where 


S (T)  - SSm 

V V 


f 

p 

f 

fll  J 

- &1Z  1 

SLII 

rii 

bJi 

yi 

(b2X2)  2 

q2 

lcfJ 

X1 

hz\ 

al2 

t 

h2Xl 


f,  II, 
2X2(T-TSLi 


SL 


<b2X2>2 


}• 


Again,  an  upper  bound  on  xmt  is  given  by  a . 2/ (a2b^) • may  he  shown  that 


II 


,11 


S (x)  < 0 for  all  x > xCT  . Also,  on  P-.-  we  have 


u 


SL 


B3 


xl(t)  - 4LU  exp 


. f 
vb2X2 


[e 


-1)}, 


x2(x)  * x2, 


y^t)  * y 


SLII 
1 ’ 

SLII 


y2(t)  ■ yj"“  exp{b2x2(x-XgJ)}, 


and 


SLII  . f,  II, 

SLII  , SLII,  II,  a12y2  , b2X2  T_TSL 

px(T)  ■ Px  exp{-a11y1  (t-xsl) j — [e 


, f 
b2X2 


-1]}  with  p 


SLII 


rx. 


SLII 


! 

! 

f 

p 

f 

n 

1 , , SLII  , SLII 

P2(T)  " p2  ‘ b2y2 

[V 

|(  al2 

SLII 

[yll 

, al2 

kJ 

<b2x2^ 2 

q2 

xi 

V £ 
b2V 

, ll\ 
(t-TSL> 


. £f  ll\ 

al2  b2X2  T_TSL  | 
e ' 


(b2x2)2 


. , SLII 
’i<T)  * ’i  - an 


f 

rx, 


'•y 


q2(t) 


m- 


12 


b2X2 


12  , SLII 


>2*2 


f + ^ 


>‘r 

{ 


. f,  H, 

j-b2X2(,'tSL)  1 _ 
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Table  I.  (cont.)  - 9 


i2-  'iv{ 


u*(t)  « 0 

V*(T)  - o 


for  T £ T 

V 


It  may  be  shown  that  S (x)  <0  and  S (x)  < 0 for  all  t > t . 

* u V V 


on  Pj^  we  have 


xx(t ) » xf  «xp{a1]Ly^*(t-T*)}, 


Xx^*)2  - y^*  coth(-A(x-T*)+B)  for  y (x^*)2  > a^*, 


x2(t)  x^/fl  “ ~2  x2*^T”Tv^ 


- 2 , V*.  2 V* 

for  y (x2  > = a2y2  * 


^/lT  y2*  " (x2*)2  tan(c(T-T*)+D)  for  y (x^V  < a^*, 


yi< T)  - yf  ’ 


C#->r) 


★ 2 v*  v* 

/sinh2(-A(x-x*)+B)  for  y (x^  )2  > , 


y,(x)«<Sr/{l-^xf(x-x;)>2 


. 2 , v*x 2 v* 

for  y ^x2  > = a2y2  * 


^ <y2*  ~ 2T"  (x2*^  )/cos2(c<t_t*)+d)  for  “f-  (x2*)2  < a2y2*’ 


^2  / v*  o ^a2  v* 
where  A - y/  (x?  )2  - y—  y2  , 


Also, 
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Table  I.  (cont.)  - 10 


12.  P^:  (concluded) 


and 


"a2P2  " b2X2q2  with  q2(T=Tv)  = q2  * 

We  have  not  been  able  to  develop  solutions  in  terms  of  "elementary"  functions 
to  the  equations  for  and  q^. 


36 


Table  I.  (concluded)  - 11 


I I / u*(t)  ■ 0 

13.  t and  P : { for  T t T 

1 v*(t)  « 0 


Results  are  similar  to  those  for  P„,  above  in  12. 

B4 


T / u*(t)  * 0 T 

PB5;  { . n for  TSL*T 

'v*(t)  = 0 


Results  are  similar  to  those  for  P_,  above  in  12. 

B4 


X / U*(T)  * 1 j 2 

>.  P^, : J for  i t s x 

M ' v*(t)  -0  SL 


Results  are  similar  to  those  for  P*^  above  in  4. 
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5.2.3.  Numerical  Examples. 


A computer  program  to  calculate  numerical  values  for  information  given 
in  Table  I was  written  in  FORTRAN  for  the  IBM  360  computer.^  A plot. of  the  field 
of  extremals  (see  Figures  5,  6,  and  7 below)  is  generated  by  this  program.  The 
closed-form  analytic  results  presented  in  Table  I are  used  whenever  possible. 

Approximate  numerical  solutions  to  transcendental  equations  (for  the  determination 

* 

of,  for  example,  tu»  t , etc.)  are  developed  by  the  well-known  Newton-Raphson 
method.  In  those  cases  for  which  closed-form  solutions  are  not  available  to  the 
state  and  adjoint  equations,  a standard  fourth  order  Runge-Kutta  numerical  inte- 
gration method  is  used.  A time  step,  At,  was  used  in  these  numerical  integra- 
tions which  yielded  agreement  to  the  fifth  place  to  the  right  of  the  decimal 
place  in  test  cases  in  which  the  approximate  numerical  solution  could  be  compared 
with  the  exact  solution. 

Parameter  sets  for  the  numerical  examples  given  in  this  paper  are  shown 
in  Table  II.  For  our  problem  (2)  we  may  consider  time  to  be  an  additional  state 
variable  so  that  the  state  space  is  five  dimensional,  i.e.  the  state  variables 
are  t,  x^,  X2»  y^,  and  y^  Thus,  unfortunately,  we  cannot  graphically  depict 
the  field  of  extremal  trajectories  but  must  be  satisfied  with  viewing  "cross- 
section"  plots  of  it. 

Table  II.  Parameter  Sets  Used  to  Generate  Numerical 


Results 

Shown 

in  Figures  5 , 6 , 

and  7 . 

Parameter 

Set 

iii 

a12 

*_2 

h 

ia 

f 

ii 

f 

^2 

f 

ll 

1 

0.003 

0.006 

0.01 

0.004 

0.005 

4.0 

8.0 

8.964 

2 

0.003 

0.006 

0.01 

0.004 

0.005 

4.0 

8.0 

11.597 

The  author  would  like  to  thank  Captain  Jeffrey  L.  Ellis  (l).  S.  Army)  for  doing 
this  work.  Subsequent  computational  contributions  were  made  by  Captain  Robert 
J.  Hill,  III  (U.  S.  Army).' 
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The  most  illuminating  plot  for  gaining  insight  into  the  structure  of  the 
optimal  fire  support  strategies  for  (2)  is  that  of  extremal  trajectories  in 
terms  of  versus  backwards  time,  r.  This  is  shown  for  parameter  set 

1 in  Figure  5.  The  corresponding  strategic  variable  values  for  X and  Y 
(i.e.  u*  and  v*)  along  each  extremal  are  also  given.  Other  plots  have 
been  considered,  but  they  provide  little,  if  any,  additional  insight. 

The  most  significant  features  of  the  field  of  extremals  shown  in  Figure 
5 are  the  two  U-singular  "surfaces" : there  is  one  in  x,^  - space  in 

V-phar.e  I and  one  in  £-space  in  V-phase  II.  In  each  phase,  X uses  the 
strategy  U*  » 1 above  the  singular  "surface"  and  the  strategy  u * 0 below  it. 
Similar  to  our  discussion  in  [32],  the  singular  surfaces  are  present  i:.  the 
field  of  optimal  trajectories  so  that  the  X artillery  avoids  "overkilling" 
either  Y^  or  Y^.  This  insight  is  obvious  when  one,  for  example,  considers 


Thun,  the  r»*te  of  destruction  of  Y^  per  unit  of  X artillery  decreases  over 

time  as  the  Y^  force  level  decreases  (see  [31]  and  [32]). 

Results  for  parameter  set  2 are  shown  in  Figure  6.  There  is  a void 

(see  p.  169  and  also  p.  187  of  [21])  in  the  field  of  extremals.  This  is  because 

* 

in  backwards  time  at  the  end  xv  of  the  U-singular  subarc  in  V-phase  II,  we 
would  have  uc(x  ) (as  given  in  Figure  3)  equal  to  1.054  if  the  adjoint 
variables  were  continuous  at  x . The  following  theorem  further  explains  this 

V ' 

situation. 
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THEOREM  2:  There  can  be  no  U-singular  subarc  beginning  in 

backwards  time  at  t+  with  a.-b.y,  = a10b„y„  for 

v 11  1/1  12  2 2 

blp2(x>2  + b2q2<Tv)y2  ‘ °- 
When  a U-singular  subarc  begins  ?*t  with 

a^b^y^  * ai2b2y2’  tbere  no  discontinuity  in  the 
adjoint  variables  at  t * (i.e.  a = 0 in  (37)). 


PROOF:  Immediate  by  (27)  and  (37), 


Q.E.D. 


Additionally,  Theorem  3 gives  the  extremal  transitions  in  X's  strategy  possible 

from  the  U-singular  surface  in  V-phase  II  as  we  work  backwards  from  t*.  Thus 

* * 

since  bjP2(Tv)x2  < b2^_q2^Tv^y2  ^or  Pararaeter  set  2,  a void  would  exist  in 

* 

the  field  of  extremals  if  the  adjoint  variables  were  continuous  at  t . 


THEOREM  3:  Assume  that  there  is  no  discontinuity  in  the  adjoint 

variables  at  t = tv  with  a^^b^y1  = a]_2b2y2*  101611 

I.  if  b1P2(Tv^x2  < b2^”q2^Tv^y2’  than  We  can  °nly 
have  u*(x)  = 0 for  t 6 (t  ,Tv+6)  where  6 > 0, 

II.  if  b1P2(Tv^x2  * b2^”q2^Tv^  ^y2*  then  W6  can  haVe 

( (a)  0, 

u*(t)  = \ (b)  (l-q2y2/(P2x2))  •b2/(b1+b2)t 


\ (b)  U 
1(c)  1, 


for  T k.  (t  ,t  +<S)  where  6 > 0. 
v v 


FROOF:  (a)  When  we  are  on  the  singular  surface  in  V-phase  II  at 

then  by  (22)  and  (23)  and  the  continuity  of  the  dual  variables  we  have 

S (t=t+)  = S (t=t"1")  = 0, 


" a2b2(b1+b2)P2X2y2{u*(Tv)- 


fb2  }L  V2l, 

' P7x2J}  * 


(47) 


where  S denotes 


dS 
u 

dt  * 


(b)  Considering  a Taylor  series  expansion  about  x = x+,  we  have  by 


the  above  for  x £ x+  = x*+ 

v v 


(t-t  )2  o o 

su(T)  =— r~su(^’ 


(48) 


where  t € (t  ,t). 

v 


(c)  When  u*(t)  = 0 for  t 6 (x  ,x  +6),  then 


*+ 


Su'T‘T„  > ■ -a2<b2)‘  W2  1 - ^3? 


q2>v 


22 


< o, 


* * 


so  that  3 {.  > 0 such  that  S (x)  < 0 for  all  x (x  ,x  +6  ) . Thus,  we 
1 u v*  v 1 ’ 

can  always  have  u*  = 0 as  we  work  backwards  in  V-phase  I from  the  U-singular 
subarc  in  V-phase  II. 


(d)  Now  let  b1P2,‘Tv^x2  ^ b2^~q?/Tv^y2'  By  ^26^’  the  U-singular 

control  in  V-phase  I u<.  = (l-q2y2/ (P2X2))  ’ b^/Cb^+b^  & 1.  Thus,  the 

*+  00  *+ 

U-singular  subarc  is  possible.  When  u*(Tv  ) = 1,  then  ^u(x=xv  ) £ 0 by 

(47).  When  inequality  holds,  it  follows  that  3 5^  > 0 such  that  Sy(x)  > 0 

for  all  x <;  (x  ,x  +6,).  Clearly,  we  cannot  have  u*  = 1 if  b,p~(x  )x_  < 
v*  v 1 12  v 2 

b2(-q2(T*))y2-  Q-E.D. 

The  same  analysis  as  used  in  the  proof  of  Theorem  3 applies  on  a 
U-singular  subarc  in  V-phasc  I when  v*  = 0.  As  long  as  (27)  holds,  one 
has  three  options  similar  to  those  of  part  II  of  Theorem  3. 


43 


5.2.4.  Filling  In  a Void. 

We  have  emphasized  that  H,  £,(t),  and  ^,(t)  are  continuous  functions  of 
time  except  possibly  at  manifolds  of  discontinuity  of  both  U*  and  V*  (see  Section 
4.3  above).  From  Theorem  3 it  follows  that  a void  must  exist  in  the  field  of 
extremals  when  these  functions  are  continuous  and  ^^P2^TV^X2  < b2'-q2^Tv^y2’ 

At  V moreover,  v*  changes  (as  we  progress  backwards  in  time)  from  1 to  0 
and  u*  from  ^/(b.+bp)  to  a different  value.  Thus,  we  have  a manifold  of 
discontinuity  of  both  U*  and  V*.  Moreover,  considering  results  given  above, 
it  is  readily  shown  that  u*(t)  remains  for  increasing  t (i.e.  backwards 
time)  equal  to  zero  once  it  changes  to  zero.  Then  from  Theorems  2 and  3 it 

4c—  "ft 

follows  that  for  b^PpUy  ^x2  b2^-C*2^Tv  ^^2  tbe  c*ua^  variQkles  must  be 

* I* 

discontinuous  to  fill  in  the  void,  and  we  must  have  u*(t)  = 1 for  t < t < r 

v u 

Furthermore,  considering  Figure  6 and  considerations  "in  the  large,"  the  maoi- 
fold  of  discontinuity  must  lie  on  the  V-transition  surface. 

Thus,  we  have  established  that  for  a^b^y^  = ai2b2y2  We  ^ave 


u*(tm)  = b./Cb^), 

y*(\)  = i 


u*(t  ) = i, 


( u VV 

V*(/)  = 


It  remains  to  determine  the  function  x^( x of  (29)  so  that  and 

may  be  computed,  and  the  jumps  in  H,  £,  and  subsequently  determined  (see 
(30)  through  (33)).  It  should  be  clear  that  it  is  impossible  to  explicitly 
determine  Tv  (*,>£)•  However,  by  computation  of  five  points  on  the  V-transition 
surface,  the  desired  partial  derivatives  may  be  estimated  by  using  linear 
approximations  to  the  appropriate  directional  derivatives  and  solving  a system 
of  four  linear  equations  in  four  unknowns.  For  parameter  sel  2 (as  the  reference 


case),  this  yielded  the  fo] lowing  estimates 


LI, 


0.0000, 


if - '°-295» 
*2 


= -0.0167, 


-0.0331. 


It  is,  therefore,  convenient  to  re-write  the  jump  conditions  across 
the  manifold  of  discontinuity  of  both  U*  and  V*. 


fi<V  - MV- 


MV  - P2<V  - P 8T 


M\>  - MV  - p 3^  - oaiibi- 


MV  "MV  -p  3^+0a12b2- 

where  p and  o are  related  by  (36).  In  this  case  the  jumps  (37)  and  (38) 
in  the  switching  functions  simplify  to 


S (t+) 

U V 


^(VVzlVl  3^"b2y2  ad) 

0 an(bi)2yi  + ai2(b2)2y2  + — at aTT (’  (52) 

1 1-a2y  2 3^  “ blylX2  3^J 


S (x+)  = 
v v 


9 l 

aUa2(bl)2ylx2  )£  0 


1_a2y2  3x„ 


3t  \ 

" biyiX2  3y^J 


Since  v*(t  ) = 0,  we  must  have  S (t  ) s 0 so  that  (50)  and  (53) 
v v v 

yield  that  o i0.  It  should  be  clear  that  0 = 0 if  and  only  if  H,  £,  and 
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q are  continuous  at  t . 

v 

that  we  must  have 


For 


o > 


S (t 
u 

a 


0, 

*+. 


the  condition  that 


> 0, 


u*(i*+>  = 1 

V 


yields 


(54) 


where  Su(xv  ) is  given  by  (52)  . Although  it  cannot  in  general  be  guaranteed 
that  (54)  will  always  hold  when  a void  in  the  iield  of  extremals  such  as  that 
shown  in  Figure  6 exists,  it  should  be  clear  that  it  must  if  the  problem  (2) 
is  to  have  a solution.  The  author  conjectures  that  this  is  true.  It  is  readily 
shown  that  when  (54)  holds,  we  have 


S 

u 


) > 0, 


r,  (t*+)  < o, 

u v ’ 


and 


*+ 

S (x  ) < 0. 

V V 


(55) 


The  appropriate  value  for  o is  determined  by  "considerations  in  the 

large:"  the  structure  of  the  entire  field  of  extremals  determines  the  value  of 

I* 

this  parameter.  In  Figure  7,  we  let  x^  denote  the  backwards  time  at  which 

I * V: 

the  U-singular  subarc  is  entered  in  V-Phase  I . Corresponding  to  xu  is  o , 

<A> 

which  yields  the  first  and  second  conditions  (18)  and  (25)  (with  u i 1)  for 

O 

<(V  J ^ *J|f 

a U-singular  subarc  with  V =0  at  t >x  . For  0 <a<a  , one  uses 

u v 

* *+  X I * 

u (t)  = 1 for  t < x < x and  then  u*(t)  =0  for  t>t  . For  a>o  , the 
v u u 

U-switching  function  Su(x)  never  changes  sign  so  that  u (x)  = 1 for  all 
* 

x>xv«  Thus,  by  manipulation  of  o,  one  may  fill  in  the  void  in  the  field  of 
extremals  in  V-Phase  I.  The  resulting  field  of  extremals  is  shown  in  Figure  7. 


5.2.5.  The  Case  of  Negligible  Small  Arms  Effectiveness . 

It  seems  appropriate  to  consider  what  happens  to  the  solution  to  the 
problem  at  hand  as  the  (relative)  effectiveness  of  (small  arms)  fire  becomes 

negligible,  i.e.  as  a^  -*■  0.  Let  us  consider  (either)  Figure  5 (or  Figure  7). 
The  ll-singular  "surface"  in  V-Phase  II  has  equation  y /yv  - 2^2^al  1^1^  ’ 


Af, 


Figure  7 


Backwards  Time,  t (minutes) 


. Filled-In  Void  in  Field  of  Extremals 


For  Parameter  Set  2. 


Thus,  as  0 with  the  other  parameters  being  held  constant,  this  singular 

"surface"  appears  higher  and  higher  on  the  y^/y_  axis  in  Figure  5.  In  t he 
limit,  the  singular  surface  does  not  appear  in  the  finite  part  of  the  plane. 

Thus,  we  have  shown  that  an  optimal  strategy  in  which  a side  divides  the  fire 
of  its  supporting  weapon  system  between  the  enemy’s  primary  (infantry)  and 
supporting  systems  can  only  occur  when  the  enemy’s  infantry  has  some  fire  effec- 
tiveness (in  the  sense  of  a non-zero  Lanchester  attrition-rate  coefficient) 
against  his  infantry. 

6.  Discussion. 

In  this  paper  we  have  examined  whe  dependence  of  optimal  time-sequential 
fire-support  strategies  on  the  form  of  the  combat  attrition  model  by  considering 
a differential  game  (see  equations  (2))  with  slightly  different  combat  dynamics 
than  those  in  the  fire-support  differential  game  considered  by  Kawara  [22]  (see 
equations  (1)).  For  this  fire-support  differential  game  (2)  we  developed  first 
order  necessary  conditions  of  optimality  and  constructed  "cross-section"  pictures 
of  the  field  of  extremals.  By  comparing  and  contrasting  the  structure  of  optimal 
fire-support  strategies  for  our  problem  (2)  with  that  for  Kawara's  fire-support 
differential  game  (1),  one  begins  to  understand  the  nature  of  the  dependence  of 
optimal  strategies  on  the  combat  dynamics  by  also  comparing  and  contrasting  the 
combat  attrition  equations  fo  these  two  differential  games. 

Our  fire-support  differential  game  (2)  was  similar  to  Kawara's  problem  (1) 
(see  [22])  except  that  we  let  the  attacker's  (i.e.  X's)  attillery  produce 
"linear-law"  attrition^  against  both  the  defender's  artillery  and  also  his  infantry 
and  let  the  defender's  infantry  produce  "linear-law"  attrition  against  the 
attacker's  infantry.  As  contrasted  with  the  optimal  time-sequen-ial  fire-support 

'For  convenience  we  use  the  term  "linear-law"  attrition  to  th note  an  attrition 
process  in  which  a target -type  undergoes  attrition  at  a rate  proportional  to  the 
product  of  the  numbers  of  fircrs  and  targets  (see  [31],  [32]). 


strategies  for  Kawara's  problem  (1)  of  always  concentrating  Jl  artillery  fire 
on  first  enemy  artillery  and  then  later  enemy  infantry  (the  timing  of  the  switch 
being  force-level  independent),  for  our  problem  (2)  the  optimal  strategy  for  one 
combatant  (the  attacker,  X)  depends  directly  on  the  enemy's  force  levels  and 
is  no  longer  to  always  concentrate  all  fire  on  either  the  enemy's  primary  or 
secondary  weapon  system.  The  latter  result,  moreover,  was  shown  to  depend  on 
the  defender's  infantry  having  some  fire  effectiveness  (in  the  sense  of  a non- 
zero Lanchester  attrition-rate  coefficient)  against  the  attacker's  infantry. 

The  solution  to  (2)  is  characterized  by  the  presence  of  singular  surfaces 
(in  Issacs'  terminology  (see  (21]),  universal  surfaces  (see  also  [18])),  a differ- 
ent one  for  each  V-phase  of  battle.  When  the  battle  state  reaches  one  of  these 
surfaces,  X follows  an  optimal  strategy  of  dividing  his  artillery  fire  bettyeen 
enemy  infantry  and  artillery  in  order  to  avoid  "overkill."  Another  characteristic 
of  the  optimal  fire-support  strategies  (not  present  for  Kawara's  [22]  problem  (i;) 
is  that  X's  optimal  strategy  may  sometimes  depend  on  V'»  distribution  of 
supporting  fires.  This  behavior  occurs  on  the  singular  surfaces.  In  fact,  X 
sometimes  must  react  instantaneously  to  changes  in  Y's  fire  distribution. 

The  development  of  even  a partial  solution  to  (2)  has  involved  a solution 
phenomenon  not  previously  reported  for  Lanrhester-type  differential  games:  the 

■J*  ' 

adjoint  (or  dual)  variables  are  discontinuous  across  a manifold  of  discontinuity 
* * 

of  both  U and  v . This  manifold  of  discontinuity  exists  for  a certain  range 
of  parameter  values  in  the  solution  to  the  problem  at  hand  (2).  Furthermore, 
there  is  a military  interpretation  to  this  manifold  of  discontinuity:  if  Y£ 

concentrates  fire  on  X2  and  X2  on  Y^,  then  when  Y2  changes  to  concentrating 
all  fire  on  X^,  X must  re-evaluate  the  worth  of  a Y ^ unit  because  it  now  has 


^The  reader  should  recall  t ha t these  represent  the  marginal  ”alues  of 
3V 

i.e.  p„(t)  = - — ,--r  where  V = V(l,x,v)  denotes  '.he  value  of  the 

r2  ox2( t)  ~ ^ 

game  (sec  [14] , f 21]) . 


force  types, 
differential 
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a direct  influence  on  the  payoff.  Such  a discontinuity  in  Lie  adjoint  variables 
is  unique  to  differential  games  'see  [3],  [4])  (i.e.  it  cannot  occur  for  a one- 
sided optimal  control  problem) . 

It  should  also  be  pointed  out  that  the  presence  of  singular  (i.e.  universal) 
surfaces  in  the  solution  to  (2)  is  apparently  independent  of  the  form  of  the 
criterion  functional  (here  terminal  payoff)  and  depends  only  on  the  combat  dynamics 
For  purposes  of  comparison  we  considered  the  same  payoff  as  considered  by  Kawara 
[22].  We  also  showed  that  the  singular  (i.e.  universal)  surfaces  can  only  be 
present  in  the  solution  when  the  defender's  infantry  has  a nonzero  casualty 

producing  capability  against  X^. 

The  problem  (2)  considered  in  this  paper  has  certain  similarities  to  the 
"War  of  Attrition  and  Attack:  Second  Version"  studied  by  R.  Isaacs  (see  pp.  330- 

335  of  [21]).  We  have,  however,  developed  a much  more  complete  solution  to  our 
problem  than  that  given  in  [I’l]  for  Mengel's  problem.  Although  this  problem  (2) 
possesses  some  similarities  to  the  Lanchester-type  optimal  control  problem 
studied  by  us  in  [31],  its  solution  has  turned  out  to  be  much  more  complex.  Our 
developments  in  this  paper,  however,  have  been  significantly  helped  by  intuition 
gained  in  the  study  of  the  simpler,  one-sided  problem  (see  [32]  for  a further 
discussion) . 

As  a result  of  our  investigation  here,  we  hope  that  a better  understanding 
of  optimal  fire-support  strategies  has  been  developed.  As  is  always  the  case, 
however,  the  insights  gained  into  the  optimization  of  combat  dynamics  from  our 
study  of  the  differential  game  (2)  are  no  more  valid  than  the  combat  model  itself. 
Our  work  here  shows  that  .he  functional  forms  of  the  various  target-type  casualty 
rates  produced  by  the  artillery  essentially  determines  the  most  significant  aspects 
of  the  structure  of  the  optimal  fire-support  strategies.  Thus,  our  study  of  this 
optimization  problem  shows  the  importance  of  determining  the  appropriate  (Lanotiesl. 
type)  model  of  combat  dynamics. 
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